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The multiplier statistics of discrete and continuous nonconservative multiplicative cascade mod- 
els, employed to describe the energy cascade in fully developed turbulence, is investigated. It is found 
to be indistinguishable due to small-scale resummation and restoration of spatial homogeneity. Only 
the subclass of positively skewed weight distributions, associated to the cascade generator, quali- 
fies to describe observed unconditional as well as conditional multiplier distributions qualitatively; 
popular log-Poisson and log-stable weight distributions do not share this property. 
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I. INTRODUCTION 



The intermittent small-scale dynamics of fully developed turbulence at very large Reynolds numbers is believed to 
be more or less scale-invariant within the so-called incrtial range, the latter being bound by the integral length scale 
L and the dissipation length scale i)«L. Consequently, both experimental H] and theoretical |^],[| efforts largely 
concentrate on scaling exponents, which are deduced from structure functions of the velocity field. However, "scaling 
exponents are not everything" as they represent for sure the simplest, but at the same time also the most insensitive 
observables characterising multi-scale processes. Because of this insensitiveness any heroic effort to derive the observed 
multi-scaling directly from the Navier-Stokes equation QJ^] remains somewhat incomplete as some physics might have 
been lost with the applied approximations and assumptions. Also, due to the finiteness of the real-world inertial 
range and the limited statistical sampling, the experimental procedure to extract multi-scaling is not unequivocal and 
quoted values for scaling exponents beyond sixth order should not be overemphasised. For these reasons and also 
in view of the quest for universality in the up-to-now unreachable limit R\ — > oo, additional observables other than 
scaling exponents are needed. 

In order to explain this philosophy, "scaling exponents are not everything" , in more detail, we discuss heuristic 
multiplicative cascade processes and show that with a clever subclass of additional observables more can be learned 
about the intrinsic cascade dynamics than only from scaling exponents. Following Richardsons historic metaphor of 
large eddies decaying into successively smaller and smaller eddies, a multiplicative cascade process phenomenologically 
describes the turbulent redistribution of energy flux from large to small scales and is able to reproduce the observed 
scaling exponents of the surrogate energy dissipation field , which are related to those of the velocity structure 
functions by the Refined Similarity Hypothesis [Q. Originally intended to directly access the multiplicative weight 
distribution associated to the energy cascade generator distributions of so-called multipliers or break-up coefficients 
|^| have been extracted from very large Taylor-scale based Reynolds number atmospheric boundary layer and wind 
tunnel flows recently see also Ref. [[l^ll . It was found that within the upper part <C I < L of the inertial 

range these multiplier distributions do in fact become scale- independent. However, they do depend on the relative 
position of parent and daughter domain; this non-homogeneity of the breakup leaves a subsequent extraction of 
scaling exponents inconsistent []l3f and indicates the presence of correlations between successive multipliers, as has 
been confirmed by studying multiplier distributions conditioned on the value of the (scale-) previous multiplier |fO|,|ll| . 
The observed correlations appear to be in conflict with the simple multiplicative cascade models, where the cascade 
generator is assumed to be identical and independent at each breakup. As has been convincingly demonstrated in 
Refs. fl4|,[l5|] this apparent conflict can be resolved: introducing non-conservative cascade generators with a positively 
skewed weight distribution and restoring Euclidean homogeneity of the ultrametric cascade models, the multiplier 
distributions resulting from a discrete binary multiplicative cascade become a scale- independent fix-point due to small- 
scale resummation and are in good qualitative agreement with the experimentally deduced distributions, including 
all observed correlation effects. 

Originally not anticipated, these multiplier distributions represent one of those wanted clever subclasses of additional 
observables, from which more can be learned about the relevance of cascade processes in turbulence and about 
the intrinsic cascade dynamics than only from scaling exponents. The experimental analysis JlO| , pT| reveals that 
unconditional multiplier distributions are observed to be scale- independent within the range 500?7 < I < L/2. In view 
of the findings of Ref. [|l5| we call this the visible cascade range. Since due to small-scale resummation it takes about 
three binary scale steps for the multiplier distributions to converge to the scale-independent fix-point the lowest part 
of the true cascade range rf ~ 500ry/2 3 w 6O77 < I < L is not visible. A further conclusion from this line of reasoning 
would be that below 77* dissipation dynamics sets in and modifies the scale-independent cascade dynamics. Besides 
this interpretation about the scope of cascade processes in turbulence, also more details about the intrinsic cascade 
dynamics can be learned from multiplier distributions [jl4 15 1. Not every weight distribution, associated to a binary 
cascade generator and reproducing observed scaling exponents within experimental error bars, qualifies to yield the 
correct distributions of multipliers and their correlations. It has to come with a positive skewncss. Weight distributions 
of log-normal or certain asymmetric binomial type share this property, whereas, for example, a log-Poisson weight 
distribution |l(| does not. 

All those conclusions have been drawn from binary discrete multiplicative cascade processes in connection with 
base-two (A = 2) multiplier distributions. Multiplier distributions referring to other scale steps have not been looked 
at. In this respect it is of interest to find out whether the latter do again reveal a scale-independent fix-point 
behaviour with additional correlations or whether they are simply artefacts of the binary discrete cascade model 
implementation. Not only this aspect is of relevance to find answers, the impact of cascade model implementations 
with scale steps other than A' = 2, i.e. A' > 2 or A' — > 1, has to be studied for the multiplier phenomenology: 
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do multiplier distributions resulting from different model implementations differ or not? Put in other words: are 
different cascade model implementations distinguishable or indistinguishable on the level of multiplier observables? - 
Theoretically the scale step A' associated to the cascade generator is not fixed. Due to the present lack of a derivation 
of cascade models from the Navier-Stokes equation, no specific choice is favoured for A', except for personal taste. 
Only if the last question can be answered with yes, we can consider multiplier distributions as a clever subclass of 
additional observables. 

The organisation of the paper is as follows: In Sect. 11 we sketch the implementation of general multiplicative 
cascade processes and briefly summarize previous work on (A = 2)-multiplier distributions resulting from a discrete 
(A' = 2) model implementation. Sect. Ill concentrates on discrete (A' > 2) implementations of multiplicative cascade 
processes and presents multiplier distributions referring to scale steps X =/= 2. A quasicontinuous implementation 
of random multiplicative cascade models is used in Sect. IV to study respective multiplier distributions referring to 
various scale steps. So-called log-stable weight distributions |l7],|l8| are considered in Sect. V and, although leading to 
correct scaling exponents jl9 20 1, shown not to reproduce the observed multiplier distributions. Several conclusions 
are presented in Sect. VI, the most important being that it does not matter how multiplicative cascade models are 
implemented (discrete or continuous), since the multiplier observables are indistinguishable. 



II. IMPLEMENTATION OF DISCRETE AND CONTINUOUS MULTIPLICATIVE CASCADE MODELS 

A. Implementation 

The following one-dimensional implementation of multiplicative cascade processes is designed to model the energy 
cascade from large (L) to small (??*) length scales: 

The initial energy flux density field at the integral scale L is set equal to e(x, L) = 1 without any loss of generality. 
For an intermediate length scale rj* < I < L we introduce a so-called weight field 

oo 

q(x,l)= qk(l)0(x-kl,l) , (1) 

k— — oo 

where the index function 

eOM) = |l !°-n <; in (2) 
v ' ; \ (x < 0, x > I) w 

represents a subdomain of length I. Each random weight qk(l) is independently drawn from a scale-independent 
probability distribution p(q) with mean (q) — 1. The energy density field 

e{x,rf) = (jjqixji)^ e(x,L) (3) 

at the smallest scale rj* is then constructed as a multiplication of weight fields belonging to a hierarchy of different 
intermediate length scales h = L/(X') 1 , where (A') 11 = L/rf. 

Let us discuss the special case for A' = 2, so that it is sufficient to consider the fields only in the domain < x < L. 
The implementation (^) can then be understood in terms of a binary cascade generator: The energy density e(l) of a 
domain of length I is redistributed onto a left and right subdomain of length 1/2 with energy densities £l(7/2) = qh£{l) 
and £#(7/2) = q^e^l), where the random weights are drawn from the splitting function 

p(qL,qn) = v{qh)p{qR) ■ (4) 

The successive application of this cascade generator from the large scale L down to the small scale 77* leads to the 
field (^). - The splitting function (Q) is not in its most general form as it does not need to factorise. However, it 
has been noted in Ref. p5[ , that because fully developed turbulence is a 3 + 1 dimensional process, the observational 
reduction of a more-dimensional multiplicative cascade process to one dimension leads to a splitting function which 
factorises almost completely. 
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B. Observables 



The standard observable of multiplicative cascade processes are scaling exponents representing the intermittency 
corrections. They are deduced from a scale-invariant moment analysis of the synergetic field (||); a naive derivation 
gives 

(e(x, i m y) = {q(x, iyr = (A) T{v = (\r iT{ " ] w 

with scaling exponents 

riy) = W , (6) 

where the independence of the m weights q(x, l±), q(x,l m ) has been exploited. This deduction is not in full 
accordance with the analysis of experimentally measured fields; those are recorded at the finest resolution scale 
77 < rj* , then averaged over larger length scales I, 



e{x,l) = j / :(./. n)<l.r' . (7) 



rx+l/2 
c-l/2 

and are not identical to the model field evolved from scales L to I [Jl5|,[l7), i.e. 

s(x, I) ^ e(x, I) . (8) 

The scale dependence of moments (e(x, l) v ) of the backward field are different from (|^) and do not show rigorous scaling 
at the very small scales Only in the large-scale asymptotic scaling exponents can be extracted unambiguously 
and coincide with the expression (^). 

Originally proposed in Ref. [|| and applied in Refs. |9|-[l^| the method of random multipliers (or break-up coefficients) 
aims at directly accessing scaling exponents via respective distributions of the former. In its most general form 
multipliers are defined as 

where e(x' , 4) represents the averaged energy density over an offspring interval of length l/X and with its centre shifted 
by x' — x = (1 — A _1 )^A with respect to the parent interval. We distinguish three special cases for the shift parameter 
A: for A = 0, ±1/2 we call M (x, I; A, A) centred, right- and left-sided, respectively. - For very large Reynolds number 
turbulent flows it has been found experimentally [9|-[T]]] that the unconditional multiplier distributions p(M(l; A, A)), 
which are sampled over x, do not depend on the scale I within the regime 77 <g; I < L/2. However, a dependence on the 
shift parameter A was noticed; this "nonhomogeneity of the breakup" obscures the relationship between multiplier 
and weight distribution and the once intended extraction of scaling exponents Jl3| . Actually the A-dependence hints 
to the existence of correlations between successive multipliers. These correlations have been observed in conditional 
multiplier distributions p(M(x', l/X; A, A)\M(x, I; A, A)) with left/right-sided (A = ±1/2) @ and centred (A = 0) 
fllf shift parameter, respectively, and with the emphasis on base-two (A = 2) scale steps. 

These latest findings appear to contradict the simple multiplicative cascade processes with independent breakups. 
However, this statement is not true! The studies in Refs. [§~4|,[l5) with nonconservative binary (A' = 2) random 
multiplicative cascade models have given a new interpretation to multiplier distributions: unconditional multiplier 
distributions represent a scale-independent fix-point due to small-scale resummation; the "nonhomogeneity of the 
breakup" is also naturally explained along this line. Once the sampling of multipliers is allowed over all and not only 
dyadic positions the correct multiplier correlations have been obtained from positively skewed weight distributions; 
for example, a log-normal weight distribution does qualify whereas a log-Poisson distribution does not qualify. These 
results show that multiplier distributions do reveal more information about the underlying weight distribution than 
might be extracted from approximate scaling exponents. 
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III. MULTIPLIER STATISTICS OF DISCRETE MULTIPLICATIVE CASCADE MODELS 



A. Discrete multiplicative cascade models with binary splittings 



For discrete multiplicative cascade models with binary splittings the scale step entering into the construction of the 
energy density field (|j) is set equal to A' = 2. Furthermore, we choose n — 10. The length of the field configuration is 
set equal to < x < 10 5 L, mimicking a long time series. As weight distribution we choose the binomial distribution 

P(q) = -^—8 (1 - (1 - a x )) + (1 - (1 + a 2 )) (10) 

with parameters ot\ = 0.3, a 2 = 0.65; this choice has the correct skewness S3 = ((q — (oV) 3 ) /(((<Z — (l)) 2 )) 3 ^ 2 = 
(a 2 — ol\) I ot\ot2 > in order to reproduce conditional base-two multiplier distributions |14[| . Other, non-binomial 
weight distributions, such as of log-normal or Pearson-III type, are also able to provide a positive skewness and lead 
to more or less identical base- two multiplier distributions [[LSI , all being in qualitative accordance with experimental 
findings 1 10 ll]] . Also multiplier distributions different from base- two do not distinguish the aforementioned weight 



distributions. Hence, we will only show results for the weight distribution (|Tc|). 

Due to small-scale resummation |15j the base-two multiplier distributions are scale-invariant in the regime 877* < 
I < L/2; for 77* < I < 877* the distributions converge to the scale-independent fix-point and for scales very close to 
the integral length scale deviations from scale-invariance set in because of the finiteness of L. Figure l(a2) shows 
simulation results for left-sided multiplier distributions (A = —0.5) at scale I — L/8: the unconditional multiplier 
distribution p(M(l; 2,-0.5)) comes close to a Beta-distribution parametrisation 

rtM) = wir M "" 1(1 - M)& " 

with Pi = /?2 = 3.2, which has been extracted from very large Reynolds number atmospheric boundary layer 
turbulence pToj ] . Once conditioned on the respective parent multiplier, correlations are observed: the distribution 
p(M(l; 2, — 0.5)|0.5 < M(2l; 2,-0.5) < 1) with a large parent multiplier becomes broader than the unconditional 
distribution and its average is shifted to a value larger than 1/2. The opposite holds for the conditional distribution 
p(M(l; 2, — 0.5) |0 < M(2/;2,— 0.5) < 0.5) with a small parent multiplier; it is more narrow than the unconditional 
distribution and its average is shifted to a value smaller than 1/2. The positive correlation between left-sided daugh- 
ter and left-sided parent multiplier is a consequence of the homogeneous sampling of multipliers over all x-positions, 
which is not restricted to dyadic cascade positions. Respective multiplier distributions of right/right combinations 
are identical to these left/left combinations. A negative correlation between daughter and parent multipliers is re- 
vealed for identical left/right and right/left combinations; their respective distributions are mirror images of those 
for left/left and right/right combinations, i.e. p(M(l; 2, 0.5)|M < M(2Z;2,-0.5) <M) = p(l - M (I; 2, -0.5)|M < 
M(2l; 2, —0.5) < M). This explains that conditional distributions, which are averaged over all four daughter/parent 
combinations, are again symmetric around M(A = 2) = 1/2; these distributions have been shown in Refs. |Io|JT^ , [l5| . 

Left-sided multiplier distributions for scale steps different from A = 2 are shown in Figs. I(al+a3). For A = 4 
the unconditional multiplier distribution converges to a fix-point for 77* <C I < L/2; its maximum is shifted towards 
smaller values since its average amounts to (M(X = 4)) = 1/A = 0.25. A conditioning on a left-sided parent multiplier 
again reveals correlations between successive multipliers: the distribution p(M (I; 4, — 0.5)|0.25 < M(4l; 4, —0.5) < 1) 
is broader than the unconditional one, whereas p(MQ; 4, — 0.5)|0 < M(4Z;4, — 0.5) < 0.25) is more narrow. - Also 
multiplier distributions of base step close to one, e.g. A = 2 1 / 4 , show fix-point behaviour in the regime 77* <C I < L/2. In 
view of the discrete (A' = 2)-steps used in the cascade evolution this result appears unexpected on first sight, but again 
it is small-scale resummation, which explains this result on second thought. The unconditional left-sided distribution 
comes with an average of (M(A = 2 1 / 4 )) = 1/A = 0.84 close to one. The respective conditional distributions 
p(M(Z; 2 1 /*-0.5)|0 < M(2 1 /% 2^ _ .5) < 0.84) and p(M(l; 2 1 / 4 , -0.5)|0.84 < M(2 1 / i l; 2 1 / 4 , -0.5) < 1) both 
deviate from the unconditional distribution, but contrary to the A = 2 and A = 4 cases, the former has now become 
broader whereas the latter turns out to be more narrow. Closer inspection reveals that this tendency has been reversed 
at around A « 1.75. 

Centred multiplier distributions for scale steps A = 4, 2 and 2 1 / 4 do qualitatively show the same behaviour as the 
respective left-sided multiplier distributions depicted in Figs. I(al)-(a3). All centred distributions are more narrow 
than their respective left-sided counterparts, reflecting again the "nonhomogeneity of the breakup". This is also 
illuminated in Fig. 2, where the A-dependence of the two exponents /3i and /?2 of the Beta-distribution (|ll|), fitted 
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to the respective unconditional left-sided and centred multiplier distributions, are shown. Since according to ([LlJ) we 
have (M) = fix/ [fix + fiz), the two exponents are related by setting (M) = A -1 . 

Conditional multiplier distributions of base step different from A = 2 have not been looked at in Refs. 
Hence, we consider our respective simulation results as predictions. 



B. Discrete multiplicative cascade models with triple and quadruple splittings 

Most often discrete multiplicative cascade models are chosen with scale steps A' = 2. Now we will consider A' = 3 
and A' = 4 and investigate the resulting multiplier distributions. Within the implementation (|l|)-(0) the parameters 
are chosen as n — 7 and < x < 10 5 L. Again for demonstration, the binomial weight distribution (Toft is picked with 
new parameters ax — 0.37, «2 = 0.8 for A' = 3 and ax = 0.42, a2 = 0.92 for A' = 4. According to the expression (||) 
the two parameter settings yield an acceptable intermittency exponent r(2) = 0.236, but actually have been fitted 
to reproduce the experimentally deduced unconditional left-sided base-two (A = 2) multiplier distribution with 
Pi = fii = 3.2 in the fix-point scale regime rf -C / < L/2. For both cases the resulting unconditional and conditional 
left-sided fix-point multiplier distributions of bases A = 4, A = 2 and A = 2 1 / 4 are almost indistinguishable from the 
respective distributions shown in Figs. I(al)-(a3) for the binary discrete multiplicative cascade implementation with 
A' = 2. The same holds for centred multiplier distributions. 

For parameter choices ax, a%, resulting in a negative skewness of the binomial weight distribution, the shape 
of the resulting unconditional left-sided fix-point multiplier distribution turns out to deviate from the symmetric 
Beta-distribution parametrisation to some small, but noticeable extend. Moreover, the correct effects observed in 
the conditional multiplier distributions can not be reproduced. - On the contrary, well adjusted log-normal weight 
distributions come again with the correct positive skewness and almost identically reproduce all unconditional and 
conditional multiplier distributions derived from the good binomial weight parametrisations. 



IV. MULTIPLIER STATISTICS OF CONTINUOUS MULTIPLICATIVE CASCADE MODELS 



A. Continuous multiplicative cascade models 



A discrete implementation of random multiplicative cascade processes is often opposed by the question "why 
discrete?" . It is true that a continuous implementation leads to nicer mathematics in terms of log-infinite divisible 



and log-stable distributions [117 20,22j,E3j, but since the relationship of the random multiplicative cascade processes 



to the Navier-Stokes equation is unclear for the moment we prefer to consider discrete and continuous versions on 
an equal footing. - As far as the numerical implementation is concerned a continuous random multiplicative cascade 
process will always be quasi-continuous. Hence, Eqs. do apply with a quasi-continuous scale step A' close to 

one; we choose A' = 2 1//8 , which compared to a binary cascade corresponds to an eightfold scale-densification. Other 
parameters are set equal to n = 80 and < x < 10 5 L. Due to scale-densification the binomial weight distribution (|Io| ) 
now comes with rescaled parameters. The setting ax = 0.1, «2 = 0.22 comes from a fit of the resulting unconditional 
left-sided base-two multiplier distribution to the expression ( |TT|) with 0x = 02 = 3.2; see Fig. I(b2). Using Eq. (^) it 
also yields an intermittency exponent of r(2) = 0.25. 

Because of small-scale resummation all multiplier distributions become again scale-independent fix-points in the 
range rj* -C / < L/2; to be more concrete, the simulations reveal a lower bound of approximately 8?7*. The uncon- 
ditional left-sided multiplier distributions of base-scale A = 4, 2 and 2 1 / 4 , which are shown in Figs. I(bl)-(b3), are 
hard to distinguish from those obtained from a scale-discrete cascade implementation, which are illustrated in Figs. 
I(al)-(a3). A similar statement can be made not only about the respective left-sided conditional distributions, but 
also for all unconditional and conditional centred multiplier distributions. This finding demonstrates that multiplier 
distributions in the fix-point regime do not distinguish discrete and continuous versions of random multiplicative 
cascade processes; again, the reason for this is small-scale resummation. 

Other suitable parameter settings of the binomial weight distribution quite easily reproduce all unconditional 
multiplier distributions shown so far. Even with a different-signed skewness, as for example with a\ = 0.195, a2 = 0.1, 
the correlation effects observed in the resulting conditional multiplier distributions go in the right direction, but 
are reduced by a factor of about 2 in magnitude. - Another interesting parameter setting is «i = 1 — 71 (A') 72 , 
«2 = (A') 72 — 1. In the limit A' — > 1 + with 71 = 72 = 2/3 the weight distribution comes with a large negative skewness 
S3 = —l/y/2{\' — 1) and leads to the infinite divisible log-Poisson cascade model of Refs. All unconditional 

multiplier distributions shown so far are again reproduced, but only very weak correlation effects are observed in 
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the resulting conditional multiplier distributions. This demonstrates that although the log-Poisson cascade model is 
capable to reproduce scaling exponents with no doubt, it fails to describe the correlation systematics of multiplier 
distributions. 

These findings suggest that in order to reproduce unconditional as well as correct conditional multiplier distributions 
the weight distribution has to possess a positive skewness, irrespective of the chosen scale-step- A' implementation of the 
underlying multiplicative cascade process. For the binomial weight distribution (|l0|) we note that for the qualitatively 
best parameter choices, i.e. (a X) a 2 ) = (0.1, 0.22), (0.3, 0.65), (0.37, 0.8), (0.42, 0.92) for A' = 2 1 / 8 , 2, 3, 4, respectively, 
the resulting skewness almost coincides at a value around S3 ~ 0.8. Here a more quantitative analysis, driven by very 
large Reynolds number data, would certainly deserve future consideration. - Multiplier distributions resulting from 
a meticulously tuned log-normal weight distribution 

"^^^(-^^y) 2 ) (12) 

are also practically indistinguishable from the results shown in Fig. 1. Qualitative best parameters for the imple- 
mentations with A' = 2 1 / 8 , 2, 3, 4 are a — 0.15, 0.42, 0.55, 0.60, respectively, reflecting an approximate, but obvious 
log-normal dependence a ~ ^/log 2 A'. Contrary to the case of the binomial weight distribution, the skewness of the 
respective log- normal distributions is not identical; it increases with increasing A' and stays positive for all four cases. 



B. Continuous-turned-discrete multiplicative cascade models 



The multiplier statistics of continuously implemented multiplicative cascade processes appears to be indistinguish- 
able from their discrete counterparts. This poses the question to what degree a continuous multiplicative cascade 
process might be described by an effective discrete multiplicative cascade process. 

For the translation of a continuous into a discrete multiplicative cascade process it is important to deal with the 
forward field. In analogy to Eq. (0) the energy density field is evolved from integral down to the intermediate binary 
target length scales I and 1/2, respectively, i.e. 



e ( x ,l) = Y[q(x,k) s(x,L) 



vi=l 
f ni+d 



\i— ni +1 



with 2rf < I < L, (A')™ ! = L/l and the scale-densification d = ln2/lnA'. We restrict the x-regime to < x < 10 5 L. 
The averaging 



1 f 



x+l 

' ft Jj 



e(x' , l)dx' 



^ 2 rX+l/2 

2^ = 1 



£l(^o) = 7/ e(x',l/2)dx' , (14) 



£r(x^) = j [ + e(x',l/2)dx' 

1 Jx+l/2 



with discrete positions x = kl, < k < 10 5 L/l is principally different from the backward averaging of Eq. (Q) and 
leads to the binary weights 



<1L 



e(x,l) 



qR = ImT ' (15) 

Sampling over all discrete positions leads to a binary splitting function p c ff (q^, qji), which still might depend on the 
scale I. 
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For the numerical simulations of a concrete quasi-continuous cascade process again an eightfold scale-densification 
(d = 8) has been chosen together with the binomial weight distribution (10) with parameters a± — 0.1, o?2 = 0.22. The 
effective binary splitting function p e g(qL,q R ) with effective weights (|15|) is illustrated in Fig. 3a for an intermediate 
length scale. For rj* < I < L/2 it is practically scale-independent; only at the very large scales a weak scale-dependence 
is observed. The projected weight distributions p c ff(9i) = / Pcs(qL, qR)dq R and p c ff(<7i?) = / Pcft{qL, qn)dqL are 
identical and come close to a log-normal distribution with a positive skewness. Observe in Fig. 3b that the projected 
weight distribution is not simply a d-fold convolution of the binomial weight distribution (|l0|); the modification is a 
result of the multivariate implementation ([I]) of the weight fields. Note also, that the effective splitting function does 
not strictly factorise; a small correlation between the two effective weights and q R exists as the relative entropy 
S rc i = -/ dq L dq R p cS (q Ll q R )\n(p cS (q Ll q R )/p eS (q L )p cS (q R )) » -0.1 < turns out to be negative. 

So far the continously implemented multiplicative cascade process has been directly translated into a (A' = 2)- 
discrete multiplicative cascade process with an effective splitting function p e s(qL,qii)- For consistency we still have to 
check, whether this effective (A' = 2)-discrete multiplicative cascade process leads to identical results for the multiplier 
statistics. The discrete implementation with A' = 2 and p e g (qL, q R ) instead of Eq. (j|) leads to the multiplier 

distributions shown in Figs. I(cl)-(c3). Qualitatively there is no difference to the distributions illustrated in Figs. 
I(al)-(a3), (bl)-(b3) and quantitative differences remain very small. 



V. NO ROOM FOR LOG-STABILITY 



In this Section we want to demonstrate again, that "scaling exponents is not everything" . In view of renormalisation 
theory so-called log-stable distributions represent mathematically very attractive parametrisations for the weight 
distribution p(q), leading its way to the denotation "universal multifractals" ]l8| . In connection with the energy 



cascade in fully developed turbulence they have already been discussed in Refs. |19|2C[], where the expression 



r(i/)=r(2)^-| (16) 

has been derived for the scaling exponents @. With an intermittency exponent of r(2) = 0.22 and a stable index of 
a w 1.6 the experimentally deduced scaling exponents pl| are reproduced remarkably well. Focusing on multiplier 
distributions, we will now reveal limitations of log-stable weight distributions. 

Concentrating first on a discrete implementation of multiplicative cascade processes, where according to (^)-(^) 
the parameters are set as A' = 2, n = 10 and < x < 10 5 L, the associated log-stable weight distribution p(q) is 
constructed by noting that — Yaq ~ S a (a, (3, fi) is distributed according to a stable distribution p5[ |. The skewness 
parameter has to be set equal to (3 — 1 and the scale and shift parameters a and fj, are fixed by the scaling exponents 
t(1) = and r(2) = 0.22. For various stable parameters < a < 2 the respective log-stable weight distributions 
are illustrated in Fig. 4a. Depending on a these distributions possess different asymmetries with respect to the 
average (q) = 1; in fact the skewness S3 = ((q — (q))^)/((q — (g)) 2 ) 1 ' 5 , which should neither be confused with the 
skewness parameter (3 nor with the notation S a of the stable distributions, amounts to S3 = —0.43, 0.06, 0.63, 1.28 
for a = 1.25, 1.5, 1.75, 2.00. If a correct skewness is really a good criterion to reproduce unconditional as well as 
conditional multiplier distributions, then according to what we now know from the previous Sections it should be 
around S3 w 1 for A' = 2. This would limit the stable index a to a regime close to 2. 

The unconditional left-sided base-two multiplier distributions resulting from the log-stable weight distributions with 
a given index a are again scale-independent in the scale range 877* < I < L/2 and are illustrated in Fig. 4b. Since the 
log-stable weight distributions have already been fitted to the scaling exponent t(2), there is no further free parameter 
left to adjust to the observed unconditional left-sided multiplier distribution (|ll| ) with (3i — P2 — 3.2. For a > 1.5 
the agreement is acceptable. Notice especially the extra contributions at M = and M = 1 for small values of a. 
This is a consequence of the algebraic tails of the stable-distributions, which lead to excess contributions at q ss in 
the corresponding weight distributions; consult again Fig. 4a. 

The conditional multiplier distributions p(M(l; 2, ±0.5) |M < M(2/;2,±0.5) < M) are also shown in Fig. 4b, 
where over all four possible left/right-sided daughter/parent multiplier combinations has been averaged. For a — 2 
a conditioning on a small parent multiplier, i.e. (M, M) = (0,0.5), leads to a distribution more narrow than the 
unconditional distribution and a conditioning on a large parent multiplier, i.e. (M , M) = (0.5, 1), leads to a broader 
distribution. This is in qualitative agreement with the experimental findings of Ref. JTcj . For a — 1.75 this tendency 
has almost vanished, for a — 1.5 it has vanished and for a — 1.25 it has even been reversed. A very similar behaviour 
is observed in the conditional centred base-two multiplier distributions, which are depicted in Fig. 4c; only for a very 
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close to 2 they are in qualitative agreement with the experimental findings of Ref . 1 1 1 . This outcome also confirms 
our conjecture from above that a correct positive skewness is needed for the weight distribution in order to reproduce 
the multiplier distributions. 

A quasi-continuous implementation of multiplicative cascade processes with log-stable weight distributions leads 
to multiplier distributions, which are identical to those shown in Figs. 4b+c for the discrete implementation. Hence, 
we arrive at the conclusion that although scaling exponents are well fitted by log-stable weight distributions with 
index a ps 1.5 jl9|,^0) the observed systematics of the multiplier distributions rules out such a stable index value. The 
log-normal limit a ~ 2 reproduces the multiplier distributions much better. 



VI. CONCLUSIONS 



Discrete and continuous implementations of geometric multiplicative cascade processes are both able to model 
multifractal characteristics observed in the surrogate energy dissipation field of fully developed turbulence. Since 
their relation to the Navier-Stokes equation remains unclear neither form of implementation is principally favoured 
over the other. Although perhaps somewhat unexpected, this indistinguishability remains once a specific class of 
observables, going beyond scaling exponents, is employed: multiplier distributions referring to various scale-steps, 
discrete or quasi-continuous, do not care whether they result from discrete or continuous versions of nonconservative 
multiplicative cascade processes. The reason for this is that the multivariate model implementation implies a small- 
scale resummation leading to scale- independent fix-point distributions. From an analysis point of view a discrete or 
continuous multiplicative cascade process with no correlations between the various branchings appears as a continuous 
multiplicative cascade process with correlations. 

There is more to learn from these fix-point multiplier distributions: not every weight distribution associated to the 
cascade generator, which matches scaling exponents within experimental error bars, qualifies to reproduce the observed 
multiplier correlations. The weight distribution has to possess a positive skewness. For example, so-called log-stable 
distributions with a stable index not close to a — 2 do not share this property and are unable to reproduce the 
observed conditional multiplier distributions. - This leaves us with a speculation: The original idea to directly access 
scaling exponents via unconditional multiplier distributions does not work due to the observed correlations between 
multipliers. However, these correlations help to restrict the broad class of cascade weight distributions, which all 
reproduce the observed lowest-order scaling exponents, to those having a positive skewness. A further restriction of 
this subclass appears to be possible with yet additional observables to be derived from the full analytic solution of 
the multivariate cascade characteristic function pq). Then, by directly accessing the cascade weight distribution, it 
is tempting to study its dependence on the Reynolds number and on the flow configuration, and, thus, to check for 
universality. 
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(X =2 ) (X =2 1/s ) (X =2 1/s ^X'=2) 




0.5 0.5 0.5 1 

M(X) 

FIG. 1. Scale-independent left-sided (A = —0.5) multiplier distributions from a discrete (al-3), a continuous 
(bl-3), and an effective, continuous-turned-discrete (cl-3) implementation of multiplicative cascade models; see text for 
model details. Solid curve: unconditional distribution p(M(l; X, — 0.5); fully underlaid curve: conditional distribution 
p(M(l; X, -0.5)|0 < M(Xl; X, —0.5) < A -1 ); hatched curve: conditional distribution p(M(7; A, -0.5) [X' 1 < M(Xl; X, —0.5) < 1). 
The circled curve in (a2-c2) represents the parametrisation dill) with 0i = 02 = 3.2. 



11 




FIG. 2. A-dependence of the Beta-distribution- (|ll|) parameter f3i for the unconditional left-sided (open stars) and centred 
(full stars) multiplier distributions resulting from the binomial weight distribution ( |ic| ) with parameters a% = 0.3, a?2 = 0.65. 
The marked points have been calculated from least-square fits and the interpolating lines are to guide the eye. 
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FIG. 3. Effective bivariate (a) and projected (b) splitting function resulting from a continuous multiplicative cascade process 
translated into a binary discrete one. For comparison the weight distribution resulting from a naive d-fold convolution of the 
infinitesimal cascade generator is shown in (b) as appropriately rescaled bars. 
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FIG. 4. Various log-stable weight distributions (a) referring to a binary discrete multiplicative cascade process and their 
resulting left /right- sided (b) and centred (c) multiplier distributions. 
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